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We study the electronic phases of the kagome Hubbard model (KHM) in the weak coupling limit 
around van Hove filling. Through an analytic renormalization group analysis, we find that there ex- 
ists a sublattice interference mechanism where the kagome sublattice structure affects the character 
of the Fermi surface instabilities. It leads to major suppression of T c for d + id superconductivity in 
the KHM and causes an anomalous increase of T c upon addition of longer-range Hubbard interac- 
tions. We conjecture that the suppression of conventional Fermi liquid instabilities makes the KHM 
a prototype candidate for hosting exotic electronic states of matter at intermediate coupling. 
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Introduction. Understanding the variations of the crit- 
ical scale T c of unconventional, i.e. electronically medi- 
ated, superconductivity is a long-standing challenge in 
condensed matter. For the cuprates, T c is non-universal 
and has been found to depend on various quantities such 
as structural parameters [1] , number of layers [2] , Fermi 
surface topology [3], and orbital content of electrons at 
the Fermi level [I]. For the latter, the d z 2 admixture to 
the dominant d x 2_ y 2 Fermi surface character has been 
suggested as a substantial influence on T c , a motif which 
is even more visible in the iron pnictides. There, at least 
all Fe tig orbitals (d xz , d yz , and d xy ) host large portions 
of electronic states in the vicinity of the Fermi surface, 
which generically necessitates a multi-band description. 
As a consequence, universal trends of pnictide supercon- 
ductivity in terms of order parameter anisotropy and T c 
sensitively depend on the structural features which de- 
termine this orbital composition [SJ |5] . 

Multi-band descriptions are both implied due to mul- 
tiple orbitals and multiple sites associated with the unit 
cell of a given lattice. While previously mentioned su- 
perconductors are all square lattices with one single site 
per unit cell, the kagome lattice [7] possesses a minimal 
three-band model due to three sites per unit cell (inset 
Fig. [IJi). For the ka gome Hubbard model (KMH), the 
three sublattices cause fundamental problems in charac- 
terizing its preferred electronic many-body phases. In 
the strong coupling limit at half filling, the kagome spin 
model exhibits strong quantum disorder fluctuations and 
has become one of the paradigmatic models of frustrated 
magnetism [8HTT3]. While the associated Mott transi- 
tion at finite coupling might still be described within 
dynamical mean field theory [TT] , the scope of collective 
electronic phases at intermediate Hubbard strength and 
general filling is particularly challenging to investigate: 
in the same way as electronic Bloch states at the Fermi 
level can involve different orbital admixtures for the mult- 
orbital case, the electronic states in the kagome lattice 
can be differently distributed among the multiple sub- 
lattices. From a tight binding perspective (Fig.[T^i), it is 



conceivable that the filling is a sensitive parameter in the 
KHM: we find two strongly dispersive bands and one flat 
band which, for appropriate fillings, has been suggested 
to be particularly susceptible to ferromagnetism along 
Stoner's criterion [T2]. While it is an ongoing demand- 
ing effort to identify kagome lattice materials at different 
electron fillings, a promising alternative route starts to 
emerge in optical kagome lattices of ultra-cold atomic 
gases, where the optical wavelengths can be suitably ad- 
justed for fermionic isotopes such as 6 Li and 40 K [T5] . 

In this paper, we take an itinerant viewpoint on the 
KHM and study how the tight-binding kagome model 
responds to weak local and longer range Hubbard in- 
teractions. The motivation is two-fold. First, we want 
to investigate what kind of competing Fermi surface in- 
stabilities emerge in the KHM revealing the interplay of 
the sublattice structure and Fermi surface topology. We 
particularly consider the regime of the dispersive bands 
around van Hove filling, where critical scales are en- 
hanced due to large density of states and nesting becomes 
relevant (Fig. [T]). For superconductivity which we expect 
to find as the generically dominant instability channel 
for weak coupling [TJ] , the multi-dimensional irreducible 
lattice representations associated with Cq v symmetry on 
the kagome lattice which are even under inversion suggest 
the possibility of topological chiral singlet superconduct- 
ing phases [15]. Second, the weak coupling limit equips us 
with a pivotal point of the KHM parameter space which 
we can solve up to analytic precision [161 117] . This pro- 
vides a valuable starting point for subsequent effective 
studies at intermediate coupling, and as such improves 
our general understanding of the KHM. 

Main results. The sublattice structure of the kagome 
lattice has a crucial influence on the character of Fermi 
instabilities, as the sublattice distribution of electronic 
states varies along the Fermi surface (Fig. [I]). We find 
d + id superconductivity in the KHM in proximity at 
van Hove filling. This finding combined with the shape 
of the dispersive bands naively suggests a similarity to 
the honeycomb model doped to van Hove filling (Fig.pl. 
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FIG. 1. (Color online). Fermi surface properties of the kagome tight binding model at n — 5/12 total filling, (a) Band structure 
in units of t resulting from the three sublattice structure of the kagome lattice (inset). The fillings n — 3/12 and n = 5/12 
(dashed horizontal line) are located at van Hove singularities as visible in the density of state plot in (b). (c) The Fermi surface 
touches the M point of the hexagonal Brioullin zone where the DOS is maximal; its topology allows for three nesting features 
one of which is Q3 = (— ir/2, v3t/2). The colors blue, red, and green label the major sublattice occupation of the Fermi surface 
states. originate from opposite shifts of Q3 and link states of similar sublattice weights, (d) The FS labels FVI defined in 
(c) assist to read off the change of sublattice occupation weights |u s (fe)| along the Fermi surface. 



There, /-wave is preferred at fillings where there are yet 
disconnected Fermi surface whereas d + id is the leading 
instability as we find one Fermi pocket (Fig. |3^) . How- 
ever, the scales of the KHM are suppressed as compared 
to the honeycomb scenario (Fig. [3^,): the KHM exhibits a 
mechanism which we call sublattice interference affecting 
the emergence of Fermi surface instabilities, as the in- 
homogeneous sublattice distribution of Fermi level states 
causes reduced nesting effects. Furthermore, while the 
usual effect of long-range Hubbard interactions would be 
to reduce the critical scale of superconductivity [13 [18], it 
gets enhanced for the KHM as the long-range interactions 
help to relieve sublattice interference effects (Fig. [3}d). 

Model. We consider the Hamiltonian 

H = Hq + Hi n t , 

H " = t J2J2 (4s> + h - c -) + (1) 

H in t = U ^ n i,t n 'i,l + "Y X! n i,T n 3,<T'^ ( 2 ) 
* (id) 

where Cj )CT denotes the annihilation operator of an elec- 
tron at site i with spin cr, m i<7 = c\ a c irJ , and \x is the chem- 
ical potential which fixes the filling. While the local Hub- 
bard interaction of scale Uq is summed over all kagome 
sites, the sum of the tight binding model with energy 
scale t as well as the nearest neighbor interaction of scale 
U\ is summed over all neighbors on the kagome lattice 
which are on different sublattices (see inset in Fig. [l^,) . 
In order to obtain the band structure of wc perform 
the Fourier transform by dividing the kagome lattice into 
unit cells containing three sites each. This corresponds 
to quantum numbers of superlattice site i, sublattice s, 
and spin a characterizing the second quantized real space 



electron operators c ( .j a . The diagonal form of ([I]) reads 
H = £ "( k )<i,n,* c k,n,*> ( 3 ) 

where n denotes the band index and the band structure 
is shown in Fig. [T^i along with the density of states in 
Fig. [TJj. The transition from real space to momentum 
space upon Fourier transform reads 

C L- = Yl <n(*)<l,n,* expi-ikiR, + r s )), (4) 

k . 1 1 

where Ri denotes the unit cell location and r s the sub- 
lattice location within the unit cell. The core informa- 
tion which is relevant for investigating the interacting 
problem ([2| is encoded in the transformation coefficients 
u sn (k) which we call sublattice weights in the following. 
For a given band n and momentum point in the Brioullin 
zone k, the coefficients obey J2 S \ u sn(k)\ 2 = 1. 

Method. We employ perturbative renormalization 
group in the two-particle pairing channel to investigate 
the superconducting instabilities 16 . Summing over all 
diagrammatic contributions up to second order in the 
interaction scales Uq and Ui, this approach is asymp- 
totically exact for infinitesimal coupling for which the 
superconducting instabilities are found in the vicinity of 
the Fermi surface. The central quantity to compute for 
the perturbative RG is the pairing vertex 

T(k,p) = r(fe(T, —ka,pa, —pa), (5) 

for incoming particles with momentum k and —k and 
outgoing particles with momentum p and p. Due to 
spin rotational invariance, the spin index effectively drops 
out of the pairing vertex as we can always constrain 
ourselves to the S z = subblock where we obtain the 
triplet channel via antisymmetrization and the singlet 
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FIG. 2. (Color online). Fermi surface properties of the honey- 
comb tight binding model at n = 5/8. (a) Band structure in 
units of t along with the DOS (upper inset) and the two sub- 
lattice structure (lower inset). Comparing the Fermi surface 
in (b) to the scenario in Fig. [lj the Fermi surface topology 
and density state are approximately identical. Inset (b): the 
sublattice occupation along the Fermi surface is homogeneous. 



channel via symmetrization of T. For infinitesimal in- 
teractions, all relevant momenta entering the interaction 
vertex are located at the Fermi surface. From the diago- 
nal form of (|5j), we compute the superconducting insta- 
bilities and obtain the pairing vertex eigenvalues Xi |16j . 
where i = 1, 2, . . . be ordered starting by Ai as the largest 
negative eigenvalue implying the strongest SC instability 
according to T c ~ E F exp[— l/|Ai|], were E F is the Fermi 
energy. Aside from the single-particle propagators as ob- 
tained by ([3]), the key object in the diagrammatic sum- 
mation for r(fc,p) is the two-particle interaction vertex 
V(kia, k 2 a, k 3 a, k^a) with incoming particles 1 and 2 as 
well as outgoing particles 3 and 4, where k 4 is given by 
momentum conservation, and, as for T, the spin index 
will effectively be dropped in the following. 

Local Hubbard interaction. Let us consider the case 
Ui = at van Hove filling n — 5/12. The Fermi surface 
is depicted in Fig. [lj;. We only consider the interaction 
in the band at the Fermi level and hence drop the band 
index n in the following. The interaction vertex takes the 
simple form 

V(k 1 ,k 2 ,k 3 ,k 4 ) = U ^2u* s (k 1 )u* s (k 2 )u s (k 3 )u s (k 4 ). 

(6) 

From (J6|, because of the locality of Uo, the only momen- 
tum dependence is given by the sublattice distribution 
weights as defined in Q . Their evolution along the Fermi 
surface is depicted through color coding in Fig. [lj; and[lp. 
Eq. [6] looks very familiar from orbital makeup factors in 
multi-orbital systems. In our case, this role is assigned to 
the sublattice weight distribution of the kagome model. 
As in the multi-orbital case, the sublattice now affects the 
nesting enhancement of particle-hole fluctuations along 
the Fermi surface. A first guess from Fermi surface topol- 
ogy without invoking the sublattice distribution would 
suggest the nesting vectors Qi = 
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7r (1, 0) , Q3 = 7T ( — |, J . As they connect Fermi sur- 
face points with mainly different sublattice occupation, 
however, the interaction vertex ([6| will be small as it is 
diagonal in the sublattice index s. This is what we call 
sublattice interference. In fact, because of sublattice in- 
terference, the most enhanced particle hole fluctuation 
channels split into 6 different nesting vectors connecting 
equal sublattice weights. For Q 3 this corresponds to a 
shift to Q± = Q 3 ± 7T (|, j^j (Fig. lL). 

It is instructive to reconcile our findings with the Hub- 
bard model on the honeycomb lattice with two lattice 
sites per unit cell (lower inset Fig. [2]), which has been 
recently investigated via RPA, 3-patch RG, weak cou- 
pling, and functional renormalization group |15[ I16[ 1191 - 
|2"2"] . There, the tight binding band structure matches 
with the dispersive bands of the kagome lattice and al- 
lows to similarly tune the honeycomb model to the equiv- 
alent van Hove filling (Fig. [2]}) . While density of states 
(upper inset Fig. [2^,) as well as Fermi surface topology 
(Fig. [2]d) exactly match with the kagome case, the sublat- 
tice weights for the honeycomb model are homogeneous 
along the Fermi surface (inset Fig. [2]d), suggesting the 
absence of sublattice interference. 

Fig. [3^ summarizes our results for local Hubbard in- 
teractions for the kagome and honeycomb tight binding 
model. We vary the doping around van Hove filling where 
we define the valence band filling n v , i.e. the fraction 
of the partially occupied band, to enable a direct com- 
parison of both cases. The van Hove filling is located 
at n v = 1/4, the Dirac cone filling at n v = 0. For 
< n v < 1/4, the Fermi surfaces are disconnected (left 
inset Fig. [3^,) , while they form one contingent pocket for 
n v > 0.25 (right inset Fig. [3^,). For the honeycomb sce- 
nario, we find that triplet /-wave SC is preferred for the 
former (B 2 representation of Cq v symmetry group), while 
d+id-wetve is preferred for the latter (E 2 representation). 
For the kagome scenario, we only find d+id-wecve in close 
proximity to van Hove filling. In detail, for d + id we find 
two degenerate SC eigenvalues X\ 2 of d-wave symmetry, 
which then in any mean field treatment yield the prefer- 
ential topological d + id chiral superconducting state in 
order to avoid loss of condensation energy due to nodes 
which would necessarily cross with the Fermi surface [15] . 
d + id has also been obtained in variational cluster ap- 
proximation calculations |23j where, however, only local 
correlations are kept and no long-wavelength features of 
the electronic phases can be addressed. 

Aside from the suppression of /-wave in the kagome 
case, the main difference between in the kagome (fc) and 
honeycomb (h) scenario is seen in the quantitative dif- 
ference of A (Fig. [3^). At van Hove filling, A& ~ l/3Ajj. 
This illustrates at infinitesimal coupling how decisively 
sublattice interference affects the emergence of supercon- 
ductivity on the kagome lattice. 

Long-range Hubbard interactions. In the case of finite 
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FIG. 3. (Color online), (a) Critical SC scale A versus valence band filling from weak coupling for local interaction only, 
presented for the kagome scenario in Fig. [T] (blue) and the honeycomb scenario in Fig. [5] (red). Dashed line denotes /-wave, 
solid line d + id-wave. The valence band filling n v — 0.25 corresponds to a van Hove point. All scales in the kagome case are 
largely reduced as compared to the honeycomb case. Below n v , the Fermi surface consists of disconnected pieces (left inset) 
and gives sizable /-wave for the honeycomb case. Above n v > 0.25, d + id is preferred along with a drop of A for larger n v . (b) 
Relative change of A for finite Ui as compared to the Uo only case as a function of U\/Uo for both lattice scenarios at n v = 0.3. 



U\ , the interaction vertex gets significantly more compli- 
cated than for the onsite interaction scenario ([6]): mo- 
mentum dependence now originates both from the har- 
monics associated with the finite interaction range as 
well as the sublattice weights. In particular, however, 
V is not diagonal in the sublattice index anymore. We 
take a representative filling in the d + id wave regime at 
n v = 0.3 and investigate the superconducting instabili- 
ties as a function of the ratio U\/Uq (Fig. [3]d). We plot 
the ratio A/Ao where Ao is the pairing vertex eigenvalue at 
Ui = 0. As elaborated on in [17] . the generic case which 
applies to the honeycomb scenario is such that long-range 
interaction should frustrate the pairing and induce a drop 
of A, which might be tuned via the degree of external ca- 
pacitive screening of the superconducting layer [TS] . The 
KHM shows a notably different behavior, as A increases 
for longer range interactions. We can understand this 
phenomenon from the perspective of sublattice interfer- 
ence and the vertex function. As the vertex becomes 
non-diagonal in the sublattice index due to longer range 
interactions, this yields a reduction of sublattice inter- 
ference effects as particle hole fluctuations between dif- 
ferent sublattice component become sizable and allow for 
reestablishing nesting enhancement given by Fermi topol- 
ogy. Altogether, the reduction of sublattice interference 
effects overcompensates the effect of harmonic modula- 
tions due to the nearest neighbor term in ([2| and yields 
a slight enhancement of A for long range Hubbard inter- 
actions (Fig. ^jp). 

Summary and outlook. The KHM shows highly anoma- 
lous behavior in terms of weak coupling Fermi instabili- 
ties such as suppressed critical scales of superconductiv- 
ity which increase upon addition of longer range Hubbard 
interactions. While this is beyond the scope of pertur- 
bative RG at infinitesimal coupling, our findings suggest 
that the KHM will likewise exhibit anomalously reduced 



critical scales of superconducting or spin density wave in- 
stabilities at intermediate coupling. This in turn might 
show the path towards stabilizing unconventional Fermi 
surface instabilities in the kagome Hubbard model. 
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